Electron trapping by electric field reversal and Fermi mechanism 
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We investigate the existence of the electric field reversal in the negative glow of a dc discharge, 
its location, the width of the well trapping the electrons, the slow electrons scattering time, and as 
well the trapping time. Based on a stress-energy tensor analysis we show the inherent instability 
of the well. We suggest that the Fermi mechanism is a possible process for pumping out electrons 
from the through, and linking this phenomena with electrostatic plasma instabilities. A power 
law distribution function for trapped electrons is also obtained. Analytical expressions are derived 
which can be used to calculate these characteristics from geometrical dimensions and the operational 
parameters of the discharge. 
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I. INTRODUCTION 

The phenomena of field reversal of the axial electric 
field in the negative glow of a dc discharge is of great im- 
portance, since the fraction of ions returning to the cath- 
ode depends on its existence and location. Technological 
application of gas discharges, particularly to plasma dis- 
play panels, needs a better knowledge of the processes 
involved. The study of nonlocal phenomena in electron 
kinetics of collisional gas discharge plasma have shown 
that in the presence of field reversals the bulk electrons 
in the cathode plasma are clearly separated in two groups 
of slow electrons: trapped and free electrons^]. Trapped 
electrons give no contribution to the current but repre- 
sent the majority of the electron population. 

The first field reversal it was shown qualitatively to 
be located near the end of the negative glow (NG) where 
the plasma density attains the greatest magnitude. If the 
discharge length is enough, it appears a second field re- 
versal on the boundary between the Faraday dark space 
and the positive column. Also, it was shown in the pre- 
viously referred theoretical work that ions produced to 
the left of this first reversal location move to the cathode 
by ambipolar diffusion and ions generated to the right 
of this location drift to the anode. For a review see 
also 0). Those characteristic were experimentally ob- 
served by laser optogalvanic spectroscopy 

Boeuf et al 4] with a simple fluid model gave an analyt- 
ical expression of the field reversal location which showed 
to depend solely on the cathode sheath length, the gap 
length, and the ionization relaxation length. They ob- 
tained as well a simple analytical expression giving the 
fraction of ions returning to the cathode and the magni- 
tude of the plasma maximum density. 

In the present Letter we introduce a quite simple 
dielectric-like model of a plasma-sheath system. This 



approach have been addressed by other authors [3, |6( to 
explain how the electrical field inversion occurs at the 
interface between the plasma sheath and the beginning 
of the negative glow. The aim of this Letter is to ob- 
tain more information about the fundamental properties 
related to field inversion phenomena in the frame of a 
dielectric model. It is obtained a simple analytical de- 
pendence of the axial location where field reversal occurs 
in terms of macroscopic parameters. In addition, it is ob- 
tained the magnitude of the minimum electric field inside 
the through, the trapped well length, and the trapping 
time of the slow electrons into the well. We emphasize in 
particular the description of the dielectric behavior and 
do not contemplate plasma chemistry and plasma-surface 
interactions. 

The analytical results hereby obtained could be useful 
for hybrid fluid-particle models (e.g., Fiala et al. 0), 
since simple criteria can be applied to accurately remove 
electrons from the simulations. 

On the ground of the stress-energy tensor considera- 
tions it is shown the inherent instability of the field in- 
version sheath. The slow electrons distribution function 
is obtained assuming the Fermi Q mechanism responsi- 
ble for their acceleration from the trapping well. 



II. THEORETICAL MODEL 

Lets consider a plasma formed between two parallel- 
plate electrodes due to an applied dc electric field. We 
assume a planar geometry, but extension to cylindrical 
geometry is straightforward. The applied voltage is V a 
and we assume the cathode fall length is I and the nega- 
tive glow + eventually the positive column extends over 
the length Iq, such that the total length is L = I + Iq. We 
have 



-V a — IE S + IqE p , 



(1) 
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where E s and E p are, resp., the electric fields in the 
sheath and NG (possibly including the positive column). 
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At the end of the cathode sheath it must be verified the 
following boundary condition by the displacement field D 



n.(D p - D a ) = a. 



(2) 



Here, a is the surface charge density accumulated at the 
boundary surface and n is the normal to the surface. In 
more explicit form, 



£pE p — e s E s — a. 



(3) 



Here, e s and e p are, resp., the electrical permittivity of 
the sheath and the positive column. We have to solve 
the following algebraic system of equations 



IqE p + lE a = — V a , 
£pE p — e s E s = a. 

They give the electric field strength in each region 



(4) 



E s = -£(l-a- 
E P = ~\ ( I - 



i <j 

V a e s 



la 

V a e £ 



(5) 



Here, we define a = 1 — — = 



Recall that in DC 



case, E p = 1 — ^f-, and e s = Eq, with ui p denoting the 
plasma frequency and v en the electron-neutral collision 
frequency. In fact, our assumption e s — Sq is plainly jus- 
tified, since experiments have shown the occurrence of 
a significant gas heating and a corresponding gas den- 
sity reduction in the cathode fall region, mainly due to 
symmetric charge exchanges processes which lead to an 
efficient conversion of electrical energy to heavy-particle 
kinetic energy and thus to heating 

Two extreme cases can be considered: i) uj p > v en , 
implying e p < 0, meaning that t coU > T p i asma , i.e, non- 
collisional regime prevails; ii) lu p < v en , e p > 0, and then 
TcoU > Tpiasma, i-e, collisional regime dominates. 

From the above Eqs.|5]we estimate the field inversion 
should occurs for the condition 1 — l - s ^- — 0, which give 
the position on the axis where field inversion occurs: 



\o_ _ v en 



(6) 



From Eq.|B|we can resume a criteria for field reversal: it 
only occurs in the non-collisional regime; by the contrary, 
in the collisional regime and to the extent of validity of 
this simple model, no field reversal will occur, since the 
slow electrons scattering time inside the well is higher 
than the the well lifetime, and collisions (in particular, 
coulombian collisions) and trapping become competitive 
processes. A similar condition was obtained in when 
studying the effect of electron trapping in ion-wave insta- 
bility. Likewise, a self-consistent analytic model have 
shown that at at sufficiently high pressure, field reversal 
is absent. 



Due to the accumulation of slow electrons after a dis- 
tance £ c = L — l Q , real charges accumulated on a sur- 
face separating the cathode fall region from the negative 
glow. Naturally, it appears polarization charges on each 
side of this surface and a double layer is created with 
a surface charge —a[ < on the cathode side and o' 2 
on the anode side. But, a' = (P • n), P = £oXeE with 
e = £o(l + Xe), Xe denoting the dimensionless quantity 
called electric susceptibility. As the electric displacement 
is the same everywhere, we have D = T)i — D 2 . Thus, 
the residual (true) surface charge in between is given by 

a = -a' 1 +a' 2 . (7) 

After a straightforward but lengthy algebraic operation 
we obtain 



where 



and 



B 



a = l1 _ 1 + E^I±) +1 (_5p 



D 



eo(e s — e p ) 



We can verify that a must be equal to 

V a e 

a = a . 

2l 



(8) 



(9) 



(10) 



(11) 



Considering that a — e^XeE, we determine the minimum 
value of the electric field at the reversal point: 



E n 



Up V a 



(12) 



Here, Xe — £™ — 1, with e rw designating the relative 
permittivity of the plasma trapped in the well. From the 
above equation we can obtain a more practical expression 
for the electrical field at its minimum strength 



E n 



V a 



T V 



el 



Tlep ^ cw 'a /I Q\ 

"~ If, T7T- U-Jj 
flew ± ep ^0 



The magnitude of the minimum electric field depends 
on the length of the negative glow l§. This also means 
that without NG there is no place for field reversal, and 
also the bigger the length the minor the electric field. 
The length of the negative glow can be estimated by the 
free path length l of the fastest electrons possessing an 
energy equal to the cathode potential fall value eV a : 



lo 



eV a 



dw 



(NF(w)) 



(14) 



Here, w is the electrons kinetic energy and NF(w) is 
the stopping power. For example, for He, it is estimated 
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TABLE I: Data used for E/p = 100 V/cm/Torr. 
Cross sections and electron temperatures are taken 
from Siglo Data base, CPAT and Kinema Software, 



Gas 


T e (eV) 


a (10~ ib cm 2 ) 


Ar 


8 


4.0 


He 


35 


2.0 


O a 


6 


4.5 


N 2 


4 


9.0 


H 2 


8 


6.0 



ph = 0.02eK 0| (in cm.Torr units, with V a in Volt). 
We denote by n ew the density of trapped electrons and 
by T ew their respective temperature. Altogether, n ep and 
T ep are, resp., the electron density and electron temper- 
ature in the negative glow region. 

By other side, we can estimate the true surface charge 
density accumulated on the interface of the two regions 
by the expression 



n ep eAAS; 



A 



(15) 



Here, Q is the total charge over the cross sectional area 
where the current flows and A£ is the width of the po- 
tential well. 



A. Instability and width of the potential well 

From Eos. II II and 1 151 it is easily obtained the trapping 
well width 



eV a 



2mi vL 



(16) 



It is expected that the potential trough should have a 
characteristic width of t he ord er in between the electron 
Debye length (Arj e = J ' ^J-f ) an d the mean scattering 
length. Using Eq.EH in a He plasma and assuming V a — 
1 kV, l = 1 m and i/ en = 1.85 x 10 9 s" 1 (with T e = 0.03 
eV) at 1 Ton (n = 3.22 x 10 16 cm -3 ) we estimate A£ w 
2.6 x 10~ 3 cm, while the Debye length is X De = 2.4 x 10~ 3 
cm. So, our Eg. 1161 gives a good order of magnitude for 
the potential width, which is expected to be in fact of 
the same order of magnitude than the Debye length. 

Table I present the set of parameters used to obtain 
our estimations. We give in Table II the estimate of the 
minimum electric field attained inside the well. The first 
field reversal at £ c ~ I ng corresponds to the maximum 
density n ew 3> n ep 0,0. So, the assumed values for 
the ratio of electron temperatures and densities of the 
trapped electrons and electrons on the NG are typical 
estimates. 

It can be shown that there is no finite configuration 
of fields and plasma that can be in equilibrium without 
some external stress [l3| . Consequently, this trough is 
dim to be unstable and burst electrons periodically (or 



TABLE II: Minimum electric field at reversal point and well 
width. Conditions: He gas, p = 1 Torr, lo = 20 cm, V a = 1 
kV, S»l = o.l, ^ = 10. 



E m (V.cnr 1 ) 


AC (cm) 


< -2.5 


2.6 x 10- a 



TABLE III: Comparison between theoretical and experimen- 
tal cathode fall distance at p=l Torr, i5/p=100 V/cm/Torr. 
Experimental data are collected from Ref. 12| . 



Gas 


KT (cm) 


(,T P (cm) 


Ar 


7.40 


0.29 (Al) 


He 


1.32 


1.32 (Al) 


H 2 


0.80 


0.80 (Cu) 


m 


0.45 


0.31 (Al) 


Ne 


0.80 


0.64 (Al) 


O2 


0.30 


0.24 (Al) 



in a chaotic process), releasing the trapped electrons to 
the main plasma. This phenomena produces local pertur- 
bation in the ionization rate and the electric field giving 
rise to ionization waves (striations). In the next section, 
we will calculate the time of trapping with a simple Brow- 
nian model. 

From Eq. we calculate the cathode fall length for 
some gases. For this purpose we took He and H2 data 
as reference for atomic and molecular gases, resp. The 
orders of magnitude are the same, with the exception of 
Ar. Due to Ramsauer effect direct comparison is difficult. 

In Table III it is shown a comparison of the experimen- 
tal cathode fall distances to the theoretical prediction, as 
given by Eq. EH Taking into account the limitations of 
this model these estimates are well consistent with ex- 
perimental data [T^ . 



B. Lifetime of a slow electron in the potential well 

The trapped electrons most probably diffuse inside 
the well with a characteristic time much shorter than 
the lifetime of the through. Trapping can be avoided 
by Coulomb collisions [lOj or by the ion-wave instabil- 
ity, both probably one outcome of the stress energy un- 
balance as previously mentioned. We consider a simple 
Brownian motion model for the slow electrons to obtain 
the scattering time r, and the lifetime T of the well. A 
Fermi-like model will allow us to obtain the slow electron 
energy distribution function. 

Considering the slow electron jiggling within the well, 
the estimated scattering time is 



2? e 



(17) 



Here, V e is the electron diffusion coefficient at thermal 
velocities. 

The fluctuations arising in the plasma are due to the 
breaking of the well and we can estimate the amplitude 
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TABLE IV: Scattering time and trapping time in the well. The parameters are: E/N = 100 Td, T g = 300 K, V a = 1 kV and 
l = 0.1 m. 



Gas V e (an 2 .8-')° Venw(s~ 1 ) b A^(cm) r (s) T (s) 

"AT 2.52 x 10 b 8.10 x 10 y 1.34 x 10~ 3 7.10 x 10"" 3.97 x 10" b 

He 5.99 x 10 6 2.39 x 10 9 1.54 x 10 -2 3.95 x 10" 11 1.70 x 10~ 5 

N 2 6.11 x 10 5 6.15 x 10 9 2.32 x 10~ 3 8.81 x 10~ 12 1.64 x 10 -4 

C0 2 1.70 x 10 6 3.60 x 10 9 6.78 x 10~ 3 2.70 x 10~ n 5.90 x 10~ 5 



"Data obtained through resolution of the homogeneous electron 
Boltzmann equation with two term expansion of the distribution 
function in spherical harmonics, M. J. Pinheiro and J. Loureiro, J. 
Phys. D.: Appl. Phys. 35 1 (2002) 

b Same remark as in a 



of the fluctuating field by means of Eq. ED We obtain 



SE„ 
Then, we have 



2 J2^- 



SE n 

E m 



lo ' 



(18) 



(19) 



In Table IV we summarize scattering and trapping 
times for a few gases. 



C. Power-law slow electrons distribution function 

As slow electrons are trapped by the electric field inver- 
sion, some process must be at work to pull them out from 
the well. We suggest that fluctuations of the electric field 
in the plasma (with order of magnitude of £ c )act over 
electrons giving energy to the slow ones, which collide 
with those irregularities as with heavy particles. From 
this mechanism it results a gain of energy as well a loss. 
This model was first advanced by E. Fermi |8j] when de- 
veloping a theory of the origin of cosmic radiation. We 
shall focus here on the rate at which energy is acquired. 

The average energy gain per collision by the trapped 
electrons (in order of magnitude) is given by 

Aw = Uw(t), (20) 

with U = £c and where w is their kinetic energy. After 
N collisions the electrons energy will be 



:ut\ 

w(t) = e t exp [ — I 



(21) 



with e t being their thermal energy, typical of slow elec- 
trons. The time between scattering collisions is r. As- 
suming a Poisson distribution P{t) for electrons escaping 
from the trapping, then we state 



P{t) = exp(-t/r)dt/T. 



(22) 



The probability distribution of the energy gained is a 
function of one random variable (the energy), such as 



100- 







o 



o 



D 
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o 



CO 0,01 




w/eV 



FIG. 1: Slow electrons distribution function vs. energy, for 
the same conditions as presented in Table IV. Solid curve: Ar, 
broken curve: N2. 



This density f w (w) can be determined in terms of the 
density P(t). Denoting by t\ = T the real root of the 
equation w = w{t\ = T), then it can be readily shown 
that slow electrons obey in fact to the following power- 
law distribution function 



f w (w)dw 



t dw 
Uf £ t w i+t/ut ■ 



(24) 



f w (w)dw = P{w < w < w + dw}. 



(23) 



Like many man made and naturally occurring phenom- 
ena (e.g., earthquakes magnitude, distribution of in- 
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come), it is expected the trapped electron distribution 
function to be a power-law (see Eq. l24f> . hence 1 + = 

n, with n — 2 -j- 4 as a reasonable guess. Hence, we 
estimate the trapping time to be 



Fig.l shows the slow electrons distribution function 
pumped out from the well for two cases: Ar (solid curve), 
and N 2 (broken curve). It was chosen a power exponent 
n = 2. Those distributions show that the higher confin- 
ing time is associated with less slow electrons present in 
the well. When the width of the well increases (from solid 
to broken curve) the scattering time become longer, and 
as well the confining time, due to a decrease of the relative 
number of slow electrons per given energy. This mecha- 
nism of pumping out of slow (trapped) electrons from the 
well can possibly explains the generation of electrostatic 
plasma instabilities. 

Note that the trapping time is, in fact, proportional to 
the length of the NG and inversely proportional to the 
electrons diffusion coefficient at thermal energies: 

T«JL. (26) 

The survival frequency of trapped electrons is v t = l/T. 
As the electrons diffusion coefficient are typically higher 
in atomic gases, it is natural to expect plasma instabil- 
ities and waves with higher frequencies in atomic gases. 
This result is in agreement with a kinetic analysis of in- 
stabilities in microwave discharges ^4|. I n addition, the 



length of the NG will influence the magnitude of the fre- 
quencies registered by the instabilities, since wavelengths 
have more or less space to build-up. Table lTvl summarizes 
the previous results for some atomic and molecular gases. 
The transport parameters used therefor where calculated 
by solving the electron Boltzmann equation, under the 
two-term approximation, in a steady-state Townsend dis- 
charge [l5j 



III. CONCLUSION 

We have shown in the framework of a simple dielectric 
model that the magnitude of the minimum electric field 
(on the edge of the negative glow) depends directly on 
the applied voltage and is inversely proportional to the 
NG length. 

The width of the well trapping the slow electrons is 
directly dependent on the applied electric field and is in- 
versely proportional to the square of the electron-neutral 
collision frequency for slow electrons. It is, as well, in- 
versely proportional to the NG length, and has typically 
the extension of a Debye length. We state that for typical 
conditions of a low-pressure glow-discharge, field reversal 
occurs whenever uj p > v en , due to a lack of collisions nec- 
essary to pump out electrons from the well. Furthermore, 
the analytical expressions obtained for the scattering and 
trapping time of the slow electrons are potentially useful 
in hybrid fluid-particle plasma modelling. 
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